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Abstract 

A study on the occurrence of Aspergillus section Nigri species on grapes from four traditional grape-producing areas in 
Greece during the 201 1/2012 vintage, and their capability to produce OTA was conducted. One hundred and twenty-eight 
black aspergilli isolates were characterised at the species level initially by the use of morphological criteria in accordance 
with appropriate keys, followed by molecular characterisation performed with Polymerase Chain Reaction-Restriction 
Fragment Length Polymorphism (PCR-RFLP) of the 5.8 ribosomal RNA gene Internal Transcribed Spacer region (5.8 rRNA 
ITS). Restriction enzyme digestion of the ITS amplicons using the Hha\, Hinfi and Rsa\, endonucleases distinguished eleven 
different patterns of restriction fragment length polymorphism (RFLP), four for each of the Hha\ and Rsal digests and three 
for Hinft. From a total number of 128 individual isolates, 124 were classified into four Aspergillus species corresponding to A. 
carbonarius, A. tubingensis, A. japonicus and A. ibericus, and the remaining 4 were classified as members of the A. niger 
aggregate. A. carbonarius and A. tubingensis being the main representative species were equally counted, with higher 
geographical representation of the former in southern and the latter in northern regions, respectively. All isolates were 
tested for their ochratoxigenic potential by use of High Performance Liquid Chromatography (HPLC) and Enzyme Linked 
Immuno Sorbent Assay (ELISA), resulting in significant interspecies differences in OTA production. 
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Introduction 

Ochratoxin A (OTA) is a naturally occurring mycotoxin, 
produced principally by a wide range oi Aspergillus and Pmicillium 
species [1] while various nephrotoxic, carcinogenic, immunotoxic, 
genotoxic and teratogenic effects [2,3,4,5] as well as Balkan 
Endemic Nephropathy [6,7] have been attributed to this 
mycotoxin. 

OTA is receiving attention due to its high incidence in a wide 
range of food commodities such as cereal-based products, coffee, 
spices, nuts, olives and grape-derived products [8,9,10,11,12,13]. 
RASFF Annual Report places OTA as the mycotoxin with most 
notifications for fruit and vegetables in 2012, while is being always 
second after aflatoxins in notifications on mycotoxins for food and 
feed for over the last decade [14]. Several pubUcations report the 
high occurrence of OTA in wine and grape products from 
different European countries, with higher OTA levels reported for 
products originating from southern European regions and in 
particular in regions with a Mediterranean climate 
[15,16,17,18,19,20]. Since 2005 the European Commission has 
imposed regulatory limits for OTA, and has established a 2 
maximum level of OTA in wine and grape products [21,22]. 

Several studies performed worldwide have shown that OTA is 
produced during infection of grapes in vineyards mainly by 



mycotoxigenic strains of black aspergilli (section Nigri), in 
particular Aspergillus carbonarius and species belonging to the 
Aspergillus niger aggregate such as Aspergillus niger and Aspergillus 
tubingensis, [18,19,20,23,24,25]. A. carbonarius although less com- 
mon than other black aspergilli, is considered to be the 
predominant species responsible for OTA contamination in 
grapes and wine, because of the ability of almost all of its strains 
to produce high levels of the toxin [26,27]. Due to the variations in 
the toxin production potency of different aspergilli, it is of great 
importance to identify the Aspergillus species accurately in order to 
define potential toxicological risks at an early stage [28] . 

Black aspergilli are difficult to classify following only morpho- 
logical criteria [29] . Though A. carbonarius can be easily recognized, 
closely related morphospecies like the biseriate species included in 
the A. niger aggregate, or uniseriate species such as Aspergillus 
aculeatus and Aspergillus japonicus, have been always difficult to 
distinguish. Among the molecular approaches used to decipher the 
Aspergillus taxonomy [30,31,32,33,34,35,36,37,38,39], Polymerase 
Chain Reaction - Restriction Fragment Length Polymorphism 
(PCR-RFLP) has been used successfully to identify Aspergillus 
species [33,39]. The non-coding Internal Transcribed Spacer 
regions (ITS) due to their high degree of sequence polymorphism 
have been used widely for fungi molecular systematics at the 
species level [40,41,42]. Standard primer pairs have been used for 
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selective amplification of these fungal sequences and have been 
proved useful for the identification oi Aspergillus species [33,39]. 

The objective of this study was to characterize (i) the biodiversity 
of black aspergiUi isolates from vineyards of four traditional grape 
producing regions in Greece during the 2011/2012 vintage using 
molecular techniques and (ii) study their ochratoxigenic potential. 
To our knowledge this is the first study with the objective to 
monitor the population of black aspergiUi by molecular means in 
Greece. 

Materials and Methods 

Ethics Statement 

Part of this study was carried out in private lands. No specific 
permissions were required for these locations/activities that took 
place and all the land owners agreed to perform this research in 
their site. The field studies did not involve endangered or 
protected species. 

Sampling and culturing conditions 

Thirty-three vineyards, located at four traditional grape- 
producing regions of Greece, were chosen for sampling during 
the 2012 harvesting period (from August to October). The 
sampled vineyards were located in Iraklion prefecture of Crete 
(n = 8), Mesogeia province of Attica (n=ll), Corinthia and 
Arcadia prefectures of Peloponnese (n = 7), Thessaloniki, Imathia, 
Fiorina and PeUa prefectures of Macedonia (n = 7). Selection of 
sampling areas was made with provision to represent the whole 
range of the typical climatic conditions of Greece, i.e. the 
Peloponnese from Southern Greece and Attica from Central 
Greece, having both a landlocked Mediterranean climate with 
vineyards of high and low altitudes, respectively, Crete from the 
Southern Aegean Sea with insular Mediterranean climate, and 
Macedonia from Northern Greece having a typical mountainous 
climate with vineyards of high altitude. 

From every vineyard five plants were selected along two major 
diagonal transects and two bunches were collected from each 
plant, resulting in a total of 10 bunches per sampled vineyard. 
Bunches were kept separate in sterile plastic bags and stored in 
portable refrigerators during transfer to the laboratory where the 
analysis took place within 24 hours from sampling. 

From each bunch five healthy berries were randomly selected 
and placed directly on the surface of Dichloran Rose Bengal 
Chloramphenicol (DRBC) medium (LabM, UK). The plates were 
incubated in the dark at 25°C for 7 days. 

Mycoflora enumeration and black aspergiUi isolation and 
identification 

After incubation, incidence of infected with black aspergiUi 
berries and distribution of most dominant grape mycoflora at 
genus level, were recorded. A representative number from the 
black aspergUli from every vineyard were isolated from DRBC and 
sub-cultured on Malt Extract Agar (MEA, LabM) and Czapek 
Dox Agar (CD, Oxoid) for further identification at species level. 
The initial identification of the different strains oi Aspergillus section 
JVigri was performed using macroscopic and microscopic morpho- 
logical criteria in accordance with appropriate keys [29,43,44]. 
The reference strains of A. carbonarius, A. niger, A. tubingensis, A. 
ochraceus and A. westerdijkiae were kindly provided by Prof N. 
Magan from the Mycology Group of Cranfield University. The 
isolates were preserved at — 80°C in the culture coUection of the 
Department of Food Science and Human Nutrition of the 
Agricultural University of Athens, Greece. 



DNA extraction 

For DNA extraction aU strains were grown in Yeast Extract 
Sucrose broth (YES; contained per litre 20 g yeast extract and 
150 g sucrose) at 30°C for 2 days. Mycelia were collected, washed 
briefly with ethanol 96%, and dried using Whatman No. 1 filter 
paper. Approximately 200 mg of mycelia from each individual 
strain were frozen in liquid nitrogen and ground to a fine powder. 
DNA extractions were performed using the NucleoSpin Plant 
DNA kit (Macherey-Nagel, Germany) according to the manufac- 
turer's instructions. 

PGR reactions and Restriction Endonuclease DNA 
digestions 

The 5.8 S-ITS region was amplified by PCR using universal 
primers ITS 1 and ITS 4 [45] . PCR reactions were performed in a 
50 /iL final volume, containing 1 x standard reaction buffer (New 
England Biolabs, UK), 2.0 mM MgClj, 300 /(M dNTPs (each), 
300 nM primers (each), 100 ng DNA template and 1.25 U of Taq 
DNA polymerase (New England Biolabs, UK). The PCR reactions 
were performed in a MJ Research PTC-200 thermal cycler (Bio- 
Rad Laboratories, USA), starting with an initial denaturation step 
at 95°C for 5 min, followed by 37 cycles consisting of 30 sec at 
95°C, 30 sec at 52°C and 40 sec at 72°C, and a final extension 
step at 72°C for 10 min. PCR products were digested with the 
Hhal, Hinjl and foal (New England Biolabs, UK) restriction 
endonucleases. Digestions were performed at 37°C for 3 h, in a 
20 /(L reaction volume containing 2 /iL of lOX reaction buffer, 
10 /(L amplicon, 1.5 U restriction endonuclease, and 0.2 /.(L BSA 
(10 |lg jjLT ) for Hhal digestions. PCR amplicons and their 
restriction digestion fragments were separated by electrophoresis 
at 100 V, 1 xTAE buffer, in 1% and 3% agarose gels, respectively. 
Agarose gels were subsequently stained in ethidium bromide 
solution (0.5 mg ml '), and DNA bands were \'isualized under 
Ultra Violet (UV) light using a Gel Doc XR^ system (Bio-Rad 
Laboratories, USA). Molecular sizes were estimated by compar- 
ison with the DNA standard GeneRuler 50 bp and 100 bp DNA 
ladders (Thermo Scientific, USA). 

Sequencing and Phylogenetic Analysis 

PCR products were purified using a QIAquick PCR purification 
kit (Qiagen, USA). Sequencing was performed for both strands 
using primers ITSl and ITS4 with the BigDye Terminator v3.1 
cycle sequencing kit (Life Technologies, USA) in an ABI3730 xl 
Genetic Analyzer (Applied Biosystems, Life Technologies, USA) 
automatic DNA sequencer (Cemia, Greece). 

Alignment of the 5.8S-ITS region sequences was performed 
using the CLUSTALQ multiple sequence alignment program 
avaUable at EBl (EMBL, UK). The ioWoWmg Aspergillus spp. 5.8S- 
ITS sequences: FJ450778, A. acidus; KF815051, A. aculeatinus; 
AJ279997, A. aculeatus; AM087614, A. awamon; AJ280010, A. 
brasihensis; AF459734, A. carbonarius; FJ491684, A. coreanus; 
FJ629326, A. costancaensis; AY656631, A. ellipticus; KC 796400, A. 
eucalypticola; HE8 18074, A. Jijiensis; AY373850, A. foetidus; 
AJ280013, A. heteromorphus; FJ629334, A. homomorphus; AY656624, 
A. ibericus; FJ629335, A. japonicus; FJ629336, A. lacticoffeatus; 
KC796401, A. neoniger, AJ223852, A. niger, FJ629352, A. piperis; 
EU159216, A. sclerotiicarbonarius; HM853552 A. saccharolyticus; 
FJ629353 A. sclerotioniger; FJ629367, A. tubingensis; AM745757, A. 
uvarumn; FJ629368, A. vadensis and FR733805, A. violaceofuscus, were 
included as outgroups for the phylogenetic analysis. Genetic 
distances were calculated using the Jukes-Cantor, parameter 
model and the phylogenetic inference was obtained by the 
Neighbour-Joining (NJ) method [46]. The NJ tree and the 
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statistical confidence of a particular group of sequences in the tree, 
evaluated by bootstrap test (1000 pseudoreplicates), were per- 
formed using the computer program MEGA version 3.0 [47] . 

OTA extraction 

All 128 black aspergiUi isolates were centre inoculated to CYA 
medium (which contained per litre: K2HPO, 1 g; Czapek 
concentrate, 10 ml; trace metal solution, 1 ml; yeast extract, 5 g; 
sucrose, 30 g; agar, 15 g) [41] and incubated at 25°C for 7 days in 
order to assess their ochratoxigenic potential. Aspergillus section 
Migri isolates was subjected to Ocliratoxin A determination 
according to the Bragulat et al. protocol [48] with a slight 
modification in the OTA extraction step. Instead of removing 3 
small agar plugs, the whole content of the Petri dish was removed 
and used for OTA extraction [49,50]. Specifically, the content of 
the Petri dish (substrate and mycelium) was weighted in order to 
express OTA production per g of substrate and extraction took 
place with 100 ml of an 80/20 methanol/water solution of HPLC 
grade purity when prepared for HPLC analysis and a 50/50 
solution when prepared for the ELISA method. The weighted 
substrates were blended with the solutions for 2 min and left for a 
total of 30 min before filtered, at first through a Whatman No 1 
filter paper, and subsequently through Millex nylon membrane 
filter of 0.2 jim pore size (EMD MiUipore Corp. BiUerica, USA), 
and kept at — 80°C until analysis. Additionally, known concen- 
trations of OTA (50, 100 & 500 ppb) were spiked on CYA and 
recovery rates for both HPLC and ELISA methods were 
estimated, resulting in satisfactory recovery percentages of 96- 
99% and 84—95% for the former and latter method, respectively. 

HPLC Analysis 

Ocliratoxin A analysis was performed using reverse-phase High 
Performance Liquid Chromatography with fluorometric detection 
(HPLC-FLD). This consisted of a JASCO AS-2055Plus auto- 
sampler, aJASCO LC-Net II/ADC system controller, a JASCO 
PU-980/LC-980-02 pump, and a JASCO FP-2020Plus fluores- 
cence detector JASCO Inc., Easton, USA). The samples were 
separated using a CI8 analytical column (250x4.6 nm, 4 /.im, 
Resteck Co., Pinnacle II, Bellefonte, USA) under isocratic 
conditions at a flow rate of 1 ml niin~' of the mobile phase 
(water/acetonitrUe/acetic acid: 99/99/2). All chemicals used for 
HPLC analysis were HPLC grade (methanol and acetic acid: 
Sigma-Aldrich Co., Germany; water and acetonitrUe: Carlo Erba 
Reactifs SDS, Val de RuHl, France). An excitation wavelength of 
333 nm and an emission wavelength at 460 nm were used for UV 
detection. Standard solutions were made from stock ochratoxin A 
solution (10.06 /ig ml ' in acetonitrile; Biopure, Romer Labs 
Diagnostics GmbH, TuUn, Austria) in mobile phase and a 
recovery study took place by spiking known concentration 
solutions to substrate and following the same extraction procedure 
as for samples. Run time for samples was 14 min with OTA being 
detected at about 1 1 min. The limit of quantification was 2.0 ng 
OTA g~' CYA (ppb), while the limit of detection was 1.0 ng OTA 
g-' CYA. 

ELISA 

Ochratoxin A quantification was additionally performed with a 
competitive direct enzyme-linked immunosorbent assay (CD 
ELISA). Veratox ELISA quantitative test kit (Neogen Corp. Ltd, 
Lansing, USA) was applied to the samples, having a detection limit 
of 1 ppb (ng g ') and range of quantification of 2-25 ppb. ELISA 
was strictiy performed with accordance to the manufacturer's 
protocol and optical densities were determined at 650 nm 
absorbance with a spectrophotometer (Synergy HT, Biotek, 



USA). Qviantification of OTA concentration was estimated with 
the aid of standard curves obtained from standard solutions 
provided within the test kit. When the determination of OTA by 
HPLC revealed concentrations greater than the upper limit of 
quantification by the ELISA method, serial dilutions of the 
extracts with the extraction solution took place in order to 
accomplish the ELISA method without extrapolation of its 
standard curves. Finally, non-ochratoxigenic isolates were omitted 
from screening with the ELISA method. 

Results 

Diversity of grapes' fungal population and Aspergillus 
section Nigri presence during the harvesting period 

The genera of filamentous fungi most abundant on sampled 
grapes were, in descending order, Aspergillus (39%), Fusarium (15%), 
Altemaria (12%), Cladosporium (11%), Rhizopus (8%), Penicillium (7%) 
and Botrytis spp. (5)%. With respect to the different regions studied, 
Aspergillus was the most frequentiy isolated species in all regions 
(32%— 55%) with the exception of Macedonia in Northern Greece 
where Altemaria was dominant (36%). 

Among the potential OTA-produciiig fungi, only black asper- 
gilli [Aspergillus section Mign group) were isolated for further 
identification to species level. Initially, the 128 different black 
aspergiUi species isolated from grapes were classified according to 
their morphological characteristics into three subgroups namely, 
Aspergillus niger aggregate, uniseriate species and A. carhonarius, and 
thereafter, a partial molecular identification took place. 

Incidence of black aspergiUi isolation from the sampled 
vineyards is presented in Figure 1. Regarding the presence of 
Aspergillus section Migri, Macedonia presented very low numbers, 
namely 47% incidence (i.e., % of contaminated out of total 
analysed grapes) with 18% distribution in the isolated mycoflora 
(i.e., % of group in total recorded fungi), while between the other 
regions occurrence of black aspergiUi was similar and always the 
dominant genera. The highest incidence and distribution were 
observed in Attica (99% incidence/55% distribution), foUowed by 
the Peloponnese (96%/34%) and Crete (76%/32%). 

Regarding A. carhonarius, the highest distribution occurred for 
isolates from Crete, with 83% of black aspergUli identified, 
followed by the Peloponnese and Attica with 45% and 32%, 
respectively, while A. carhonarius isolation was most scarce in 
vineyards of Macedonia (15%). As regards uniseriate spp. a sole 
isolate was identified originating from Attica prefecture. According 
to the spore size that can discriminate between black uniseriate 
aspergiUi, the isolate could be characterized as A. japonicus. 

Molecular characterization Aspergillus isolates by PCR- 
RFLP 

In order to monitor the representation of black aspergiUi [A. 
aculeatus, A. awamori, A. hrasiliensis, A. carhonarius, A. foetidus, A. 
ihericus, A. japonicus, A. niger, A. tuhingensis and^. uvarumn) reported to 
be isolated from grapes [20,25,34,37,38] we focused on a partial 
molecular characterization of the 128 isolates at the species level 
by use of PCR-RFLP of the 5.8S-ITS region. PGR mediated 
amplification of the ribosomal region was performed using the 
ITSl and ITS4 primers resulting to amplicon molecular sizes of 
approximately 600 bp for aU samples. RFLP analysis was done 
using the Hhal, Hinjl and Rsal restriction endonucleases in order to 
differentiate among species. The restriction patterns of reference 
strains as presented in Figure 2 were used for comparison with the 
individual isolates restriction digestion patterns. Hhal and HinJl 
enzyme digestion can differentiate A. carhonarius from A. niger/ A. 
tuhingensis since they both generate 2 different restriction patterns 
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Figure 1. Aspergillus section Nigri strains distribution. A. section Nigri distribution in total mycoflora (bars) and different species distribution 
within the A. section Nigri group (percentages inside bars) for each sampled prefecture. The percentages underneath prefectures' names refer to 
distribution and incidence of A. section Nigri. 
doi:1 0.1 371/journal.pone.0093923.g001 



designated here as A (for A. carbonarius) and B (for A. niger and A. 
tubingensis). Rsal digestion can distinguish further A. niger from A. 
tubingensis generating different restriction patterns for these two 
species, pattern A (common for A. niger and A. carbonarius) and 
pattern B (for A. tubingensis). A total of 128 black Aspergillus isolates 
previously characterized by morphological criteria were analyzed 
by PCR-RFLP and their restriction patterns were compared with 
those obtained from the reference strains. Digestion with Hi?ifi 
resulted to 6 1 A and 66 B pattern profiles and 1 for a new profile 
designated as C and closely resembles profile B. Hhal digestion 
resulted to 60 A, 66 B pattern profiles and 2 for new pattern 
profiles designated as C and D. These new profiles for both 
enzymatic digestions can be seen for samples Aj-1 and Ai-1 in 
Figures 3 and 4. Restriction digestion with Hhal can also 
differentiate uniseriate {A. japonicus, A. aculeatus, A. uvarum) from 
all biseriate black aspergUli species [33,39]. Depending on the 
pattern which will result from the Hinjl digestion, A. japonicus can 
be further distinguished from A. aculeatus [39]. According to 
restriction endonuclease site profile study (data not shown) of 
characterized black uniseriate aspergiUi ITS sequences deposited 
in public sequence databases, Hinji digestion results in a double 
band visible pattern for either A. japonicus and A. uvarum and a three 
band visible pattern for A. aculeatus. The 66 strains which gave a B- 
B pattern combination after HinJl and Hhal digestions were 
subjected to digestion with the Rsal, resulting to 62 B and 4 A type 
restriction patterns. The strain with the C-C Hinjl-Hhal, restriction 
digestion pattern when digested with Rsal, resulted to a new profile 
designated as C. This profile closely resembles Rsal profile A and 
can be seen in Figure 4. The strain with the A-D HinJl-Hhal, 
restriction digestion pattern, when digested with Rsal resulted to a 
pattern nearly indistinguishable from profile A (Figure 4). Though 
the molecular size of the largest DNA fragment generated (slightly 
above 500 bp) runs somehow faster that the corresponding band 
from profile A, the visible difference is not so obvious. This new 
profile was designated as D. According to the PCR-RFLP results 
for the 128 isolates, 60 (46.88%) were characterized as A. 
carbonarius, 62 (48.44%) as A. tubingensis, 4 (3.12%) as members of 
the A. niger aggregate, 1 (0.78%) - resembling the pattern C-C-C 
for HinJl, Hhal and Rsal digestions - could be characterized as A. 



japonicus / A. uvarum, and 1 (0.78%) with an A-D-D Hinjl, Hhal and 
Rsal digestion pattern could not be clearly classified to any of the 
previously referred species. A detailed classification according to 
both phenotypical and molecular data can be seen in Table 1. 

Sequencing and Phylogeny study 

Representative isolates (spp. Ac-28, At-10, Ac-36, Ai-1, An-2, 
Aj-1, At-46, Ac-57, An-4 and At-59) for each species as resulted 
from PCR-RFLP analysis were subjected to sequencing for results 
verification. A part from the resulted sequences alignment can be 
seen in Figure 5. The DNA sequence restriction pattern analysis 
was in accordance with these of reference strains for all isolates 
sequenced except isolate Ai-1. Isolate Ai-1 had a Hinjl restriction 
site pattern common to A. carbonarius, however, the presence of an 
extra Hhal restriction site (nucleotides 143 to 146) generated a 
Hhal restriction site profile similar to this of A. japonicus. A 
presentation of the DNA fragments generated according to the 
Hinjl, Hhal and Rsal restriction sites, in all sequenced isolates and 
representative reference isolates of black aspergUli species isolated 
from grapes can be seen in Table 2. The fragment sizes that are 
generated resemble the sizes of the fragments visualized after gel 
electrophoresis of the amplicon digests. 

BLAST analysis of all sequences indicated that 3 clones had 
high similarity to A. carbonarius, 3 to A. tubingensis, 2 to A. niger I A. 
awamori, 1 to A. japonicus / A. aeuleatus and 1 to A. ibericus. The results 
support the characterization of the 10 isolates deduced from the 
RFLP analysis and suggest that the isolate with the A-D Hinjl-Hhal 
RFLP pattern can be classified as A. ibericus. 

The phylogenetic relationship between the isolates oi Aspergillus 
species is illustrated in an N-J method cluster analysis (Figure 6). 
The ten sequences were aligned together with reference strains 
sequences retrieved from NCBI resulting to the formation of five 
clades. Isolate Aj-1 was aligned with uniseriate A. japonicus, A. 
uvarum and A. aculeatus forming together with A. aculeatinus, A. 
Jijiensis and A. violaceqfuscus a clade separated from all other biseriate 
AspergiUi sequences with a 100% bootstrap value. Two other major 
clades were formed. The first consisted by A. carbonarius / A. ibericus 
isolates (Ac-28, Ac-36, Ac-57 and Ai-l) together with A. carbonarius, 
A. ibericus and A. bra.silien.sis outgroup reference sequences, each one 
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Figure 2. Ribosomal 5.8S-ITS region restriction digestion patterns of Aspergillus reference strains. Restriction digestion patterns 
(designated as A and B) of ribosomal 5.8S-ITS DNA amplicons from Aspergillus reference strains, after digestion with the restriction endonucleases 
Hinl\, Hha\ and Rsa\. Ac: Aspergillus carbonarius, An: Aspergillus niger, At: Aspergillus tubingensis, Aw: Aspergillus westerdijkiae, Ao: Aspergillus 
ochraceus, L: Low molecular weight DNA ladder (molecular sizes are 766, 500, 350, 300, 250, 200, 150, 100, 75, 50 and 25 bp respectively). 
doi:1 0.1 371 /journal.pone.0093923.g002 



of the last three specie.s forming a .separate .subgroup inside the 
clade. A. carbonarius sequences were aligned together and were 
separated from the A. ibericus sequences with a 98% bootstrap 
value. A. hrasilimsis formed a separate branch in this clade with a 
79% bootstrap value. The second major clade was formed by A. 
tubingensis/ A. aggregate isolates (At-10, At-46, At-59 and An-2, 
An-4 respectively) and A. niger aggregate outgroup reference 
sequences. A. tubingensis sequences were clustered together with 
other black aspergiUi (not reported to be found in grapes), forming 
a subgroup separated from all other A. niger aggregate {A. niger, A. 
awamori and A.foetidus) species of the clade, with a 86% bootstrap 
value. Two other minor clades were also formed, one by A. 



saccharoljticus together with A. homomorphus and the other by A. 
ellipticus together with A. heteroniorphus. 

Ochratoxigenic potential of Aspergillus section Nigri 
species 

A total of 128 isolates from grapes belonging to Aspergillus 
section Migri, were tested in vitro for their potential to produce 
ocliratoxin A. For the estimation of their ochratoxigenic potential 
a Competitive Direct ELISA procedure was also applied apart 
from the official HPLC method. Results showed that of the 60 A. 
carbonarius isolates 59 (98.33%) produced OTA from aU geographic 
regions assayed. In contrast, among the A. niger aggregate, only I 



SF56 AJ-1 




Hhal 



Hinfl 



Figure 3. Ribosomal 5.8S-iTS region restriction digestion patterns of Aspergillus grape isolates. Restriction digestion patterns (designated 
as A, B and C) of ribosomal 5.8S-ITS DNA amplicons from various Aspergilli grape isolates (presented as isolate designations), after digestion with the 
restriction endonucleases HhM and HinA. 50 bp (I) and 100 bp (L) DNA ladders are also shown. 
doi:1 0.1 371/journal.pone.0093923.g003 
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ABBCD ABBCA AABCD 



Hhal 



Hinn 



Rsal 



Figure 4. Restriction digestion patterns of sequenced ribosomal 5.8S-ITS region amplicons. Restriction digestion patterns (designated as 
A, B, C and D) of five sequenced ribosomal 5.8S-ITS DNA amplicons from five different Aspergilli grape isolates (presented as isolate designations), 
after digestion with the restriction endonucleases Hhal, HinA and Rsa\. Each isolate is a representative of the five different Aspergillus species 
characterized in this study. L: DNA ladder. 
doi:10.1371/journal.pone.0093923.g004 



out of 62 isolates (1.61%), identified as A. tuhingensis, had the 
potential to produce OTA. Finally, neither A. ibericus nor A. 
japonicus isolates were capable for OTA production. 

Referring to ranges of OTA production, within the A. carbonarius 
species 1 isolate (1.67%) produced less than 1 ng OTA g ' CYA, 
19 isolates (31.67%) produced between 1 and 100 ng OTA g ' 
CYA, 24 (40%) between 100 and 1000 ng OTA g"' CYA, 13 
(21.67%) between 1000 and 10000 ng OTA g"' CYA and 3 
isolates (5%) had an extremely high OTA potential with quantities 
greater than 10000 ng OTA g"' CYA (Table 3). As regards the 
single A. tubingensu OTA producing isolate it had a potential of 
227 ng OTA g"' CYA. 

Regarding the performance of the ELIS A test kit quantification, 
results were linearly correlated with those of the HPLC method 
with an r^ = 0.905 (Figure 7), but with an underestimation of the 
produced OTA at concentration levels above 5000 ppb. Results of 
both methods are presented in Table 1, together with the 
characterisation of the isolates after both phenotypic and 
molecular means. 

Discussion 

Incidence of isolation and distribution of Aspergillus section Mgri 
group among the fungal genera were in accordance with previous 
reports from Greece [20,51] and other Mediterranean countries 
[52,53,54]. A very high incidence was observed in our samples 
since 23 out of 33 vineyards sampled presented 100% presence for 
black aspergilli, although some studies have reported lower indices 
for vineyards from some European and South American countries 
[19,55,56,57]. 

Regarding the molecular characterisation of the isolated black 
aspergilli, a PCR-RFLP method was applied in order to 
dilferentiate the Aspergillus strains. RFLP analysis of the rDNA 



5.8S ITS products was performed to identify diflFerent isolates of 
black Aspergillus species from grapes compared to reference strains. 
A total of 128 black Aspergilli isolate-ITS amplicons were digested 
using the Hhal, Hinjl and Rml restriction digestion endonucleases 
resulting in 4 diflFerent RFLP patterns for Hhal and Rsal 
respectively and 3 patterns for the Hirifl enzyme digests. In 
previous studies [33,39,58,59,60] according to authors' descrip- 
tions these restriction endonucleases have been used successfully to 
differentiate the A. carbonarius, A. niger, A. tubingensis, A. japonicus and 
A. aculeatus species. It was shown that Hirtfl enzyme digestion can 
diflFerentiate A. carbonarius from A. niger/ A. tubingensis, and so does 
Hhal. Furthermore, Hhal can diflFerentiate the uniseriate A. 
japonicus, A. aculeatus and A. uvarum species from all three previously 
reported biseriate aspergilli species. Hirifl can distinguish further A. 
japonicus/A. uvarum from A. aculeatus and foal can distinguish A. 
tubingensis from all other species belonging to the A. niger aggregate 
[33,39,58,59,60]. 

According to the digestion pattern combinations resulted from 
the PCR-RFLP analysis of the 128 Aspergilli isolates and reference 
species, 60 were characterized as A. carbonarius, 62 as A. tubingensis 
and 4 as species belonging to the A. niger aggregate other than A. 
tubingensis. These results are in accordance with results presented 
by other authors, supporting the use of the Hhal, Hinjl and Rsal 
restriction endonucleases for accurate characterization of A. 
carbonarius, A. tubingensis and, A. niger species [33,39,59,60]. Two 
other digestion pattern combinations have resulted using the Hhal 
and Hinjl enzymes. One combination resumed a pattern presented 
by Spadaro et al. [39], suggesting that the isolate possibly belongs 
to the A. japonicus/A. uvarum species. The other digest pattern 
combination suggested a characterization of the isolate as A. 
carbonarius according to Hinjl digestion; however, this could not be 
supported from the Hhal digestion profile which closely resembled 
this of A. japonicus. 
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Table 1. 


Identification, origin and OTA production 


of A. section Nigri isolates. 








Spp. - No 


Species 


Origin 


OTA level (ppb) 










HPLC 


ELISA 


AM 


A. tubingensis 


Crete 


<LOD 




At-2 


A. tubingensis 


Crete 


<LOD 




Ac-1 


A. carbonarius 


Crete 


31.01 


32.20 


Ac-2 


A. carbonarius 


Crete 


27.79 


43.16 


Ac-3 


A. carbonarius 


Crete 


42.40 


44.85 


Ac-4 


A. carbonarius 


Crete 


1 208.42 


803.52 


Ac-5 


A. carbonarius 


Crete 


814.63 


665.20 


Ac-6 


A. carbonarius 


Crete 


1208.79 


906.39 


Ac-7 


A. carbonarius 


Crete 


1040.77 


550.92 


Ac-8 


A. carbonarius 


Crete 


1 290.46 


590.72 


Ac-9 


A. carbonarius 


Crete 


1270.25 


585.80 


Ac-10 


A. carbonarius 


Crete 


116.90 


62.46 


Ac-11 


A. carbonarius 


Crete 


6512.18 


4828.77 


Ac-1 2 


A. carbonarius 


Crete 


442.60 


290.79 


Ac-1 3 


A. carbonarius 


Crete 


1013.88 


405.31 


Ac-1 4 


A. carbonarius 


Crete 


1274.33 


837.26 


Ac-1 5 


A. carbonarius 


Crete 


706.68 


415.62 


Ac-1 6 


A carbonarius 


Crete 


62.44 


31.88 


Ac-1 7 


A. carbonarius 


Crete 


68.88 


13.42 


Ac-1 8 


A. carbonarius 


Crete 


347.18 


61.63 


Ac-1 9 


A. carbonarius 


Crete 


229.36 


79.67 


Ac-20 


A. carbonarius 


Crete 


397.93 


1 07.69 


Ac-21 


A. carbonarius 


Crete 


556.58 


140.77 


Ac-22 


A. carbonarius 


Crete 


44.58 


23.65 


Ac-23 


A. carbonarius 


Crete 


355.52 


279.08 


Ac-24 


A. carbonarius 


Crete 


94.86 


18.25 


Ac-25 


A. carbonarius 


Crete 


760.72 


299.86 


Ac-26 


A. carbonarius 


Crete 


644.05 


615.97 


At-3 


A. tubingensis 


Crete 


227.37 


314.39 


Ac-27 


A. carbonarius 


Crete 


<LOD 




At-4 


A. tubingensis 


Crete 


<LOD 




At-5 


A. tubingensis 


Crete 


<LOD 




At-6 


A. tubingensis 


Crete 


<LOD 




Ac-28 


A carbonarius 


Crete 


1 2782.62 


5010.56 


Ac-29 


A. carbonarius 


Crete 


34825.69 


30456.77 


Ac-30 


A. carbonarius 


Crete 


86.00 


58.42 


At-7 


A. tubingensis 


Attica 


<LOD 




Ac-31 


A. carbonarius 


Attica 


2263.18 


742.60 


At-8 


A. tubingensis 


Attica 


<LOD 




Ac-32 


A. carbonarius 


Attica 


31.68 


17.79 


At-9 


A. tubingensis 


Attica 


<LOD 




At-10 


A. tubingensis 


Attica 


<LOD 




AMI 


A. tubingensis 


Attica 


<LOD 




Ac-33 


A. carbonarius 


Attica 


18752.53 


7221.62 


Ac-34 


A. carbonarius 


Attica 


748.87 


216.94 


AM 2 


A tubingensis 


Attica 


<LOD 




AM 3 


A. tubingensis 


Attica 


<LOD 




AM 4 


A tubingensis 


Attica 


<LOD 
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Table 1. Cont. 





Spp. - No 


Species 


Origin 


OTA level (ppb) 










HPLC 


ELISA 


Ac-35 


A. carbonarius 


Attica 


553.06 


341.50 


AM 5 


A. Wbingensis 


Attica 


<LOD 




Ac-36 


A. carbonarius 


Attica 


121.32 


116.19 


An-1 


A. niger aggr. 


Attica 


<LOD 




Ac-37 


A. carbonarius 


Attica 


54.13 


69.34 


Ac-38 


A. carbonarius 


Attica 


57.92 


43.92 


AM 6 


A. Wbingensis 


Attica 


<LOD 




AM 7 


A. wbingensis 


Attica 


<LOD 




AM 8 


A. Wbingensis 


Attica 


<LOD 




AM 9 


A. Wbingensis 


Attica 


<LOD 




Ac-39 


A. carbonarius 


Attica 


618.33 


393.06 


At-20 


A. Wbingensis 


Attica 


<LOD 




At-21 


A. Wbingensis 


Attica 


<LOD 




Ac-40 


A. carbonarius 


Attica 


1 82.46 


97.34 


Ac-41 


A. carbonarius 


Attica 


1 94.28 


107.03 


Ac-42 


A. carbonarius 


Attica 


17.64 


7.85 


At-22 


A. tubingensis 


Attica 


<LOD 




At-23 


A. tubingensis 


Attica 


<LOD 




At-24 


A. tubingensis 


Attica 


<LOD 




Ac-43 


A. carbonarius 


Attica 


4.71 


7.73 


Ac-43 


A. carbonarius 


Attica 


2.07 


2.03 


At-25 


A. tubingensis 


Attica 


<LOD 




Ac-45 


A. carbonarius 


Attica 


671.44 


496.88 


AM 


A. ibericus 


Attica 


<LOD 




At-26 


A. tubingensis 


Attica 


<LOD 




At-27 


A. tubingensis 


Attica 


<LOD 




Ac-46 


A. carbonarius 


Attica 


4.26 


3.83 


An-2 


A. niger aggr. 


Attica 


<LOD 




At-28 


A. Wbingensis 


Attica 


<LOD 




At-29 


A. Wbingensis 


Attica 


<LOD 




At-30 


A. Wbingensis 


Attica 


<LOD 




Ac-47 


A. carbonarius 


Attica 


3924.93 


4041.81 


At-31 


A Wbingensis 


Attica 


<LOD 




Ac-48 


A. carbonarius 


Attica 


2725.74 


2998.55 


At-32 


A. Wbingensis 


Attica 


<LOD 




At-33 


A. Wbingensis 


Attica 


<LOD 




Ac-49 


A. carbonarius 


Attica 


23.24 


9.87 


At-34 


A. Wbingensis 


Attica 


<LOD 




At-35 


A. Wbingensis 


Attica 


<LOD 




At-36 


A. Wbingensis 


Attica 


<LOD 




At-37 


A. Wbingensis 


Attica 


<LOD 




At-38 


A. Wbingensis 


Attica 


<LOD 




At-39 


A. Wbingensis 


Attica 


<LOD 




At-40 


A. Wbingensis 


Attica 


<LOD 




AMI 


A. Wbingensis 


Attica 


<LOD 




At-42 


A. Wbingensis 


Attica 


<LOD 




Aj-1 


A. japonicus 


Attica 


<LOD 





At-43 A. Wbingensis Peloponnese <LOD 
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Table 1. Cont. 



Spp. - No 


Species 


Origin 


OTA level (ppb) 












HPLC 


ELISA 


At-44 


A. 


tubingensis 


Peloponnese 


<LOD 




At-45 


A. 


tubingensis 


Peloponnese 


<LOD 




At-46 


A. 


tubingensis 


Peloponnese 


<LOD 




At-47 


A. 


tubingensis 


Peloponnese 


<LOD 




At-48 


A. 


tubingensis 


Peloponnese 


<LOD 




Ac-50 


A. 


carbonarius 


Peloponnese 


90.34 


59.40 


Ac-51 


A. 


carbonarius 


Peloponnese 


141.51 


82.34 


Ac-52 


A. 


carbonarius 


Peloponnese 


82.00 


56.81 


Ac-53 


A. 


carbonarius 


Peloponnese 


65.44 


55.56 


An-3 


A. 


niger aggr. 


Peloponnese 


<LOD 




Ac-54 


A. 


carbonarius 


Peloponnese 


684.39 


675.38 


Ac-55 


A. 


carbonarius 


Peloponnese 


741.88 


662.32 


Ac-56 


A. 


carbonarius 


Peloponnese 


247.80 


152.97 


At-49 


A. 


tubingensis 


Peloponnese 


<LOD 




At-50 


A. 


tubingensis 


Peloponnese 


<LOD 




Ac-57 


A. 


carbonarius 


Peloponnese 


745.14 


686.79 


At-51 


A. 


tubingensis 


Peloponnese 


<LOD 




Ac-58 


A. 


carbonarius 


Peloponnese 


167.52 


116.37 


At-52 


A. 


tubingensis 


Peloponnese 


<LOD 




At-53 


A. 


tubingensis 


Macedonia 


<LOD 




At-54 


A. 


tubingensis 


Macedonia 


<LOD 




Ac-59 


A. 


carbonarius 


Macedonia 


1496.38 


605.75 


At-55 


A. 


tubingensis 


Macedonia 


<LOD 




At-56 


A. 


tubingensis 


Macedonia 


<LOD 




An-4 


A. 


niger aggr. 


Macedonia 


<LOD 




At-57 


A. 


tubingensis 


Macedonia 


<LOD 




At-58 


A. 


tubingensis 


Macedonia 


<LOD 




At-59 


A. 


tubingensis 


Macedonia 


<LOD 




At-60 


A. 


tubingensis 


Macedonia 


<LOD 




At-61 


A. 


tubingensis 


Macedonia 


<LOD 




Ac-60 


A. 


carbonarius 


Macedonia 


3654.06 


2273.10 


At-62 


A. 


tubingensis 


Macedonia 


<LOD 





Species identification and origin of A. section Nigri spp. isolated from the present study and their ochratoxigenic potential assayed by HPLC and ELISA methods. 
Limit of Quantification (LOQ) 2 ng OTA g"' CYA and Limit of Detection (LOD) 1 ng OTA g"' CYA. 
doi:l 0.1 371 /journal.pone.0093923.t001 



Sequencing of representative isolate DNAs supported the RFLP 
data since the restriction site profile of each isolate generated DNA 
fragments with molecular sizes that corresponded to the bands 
found after electrophoresis of the amplicon digests. An extra HhtA 
site is present in the sequence of isolate Ai-1 generating a 
restriction digestion profile different from the rest that are 
observed. A nucleotide substitution, C instead of T that is present 
in nearly aU other sequences generates this extra Hhal site. The 
results presented in Table 2 show that the DNA fragments 
generated according to the Hinji restriction sites are nearly equal 
in size and close to the DNA fragment sizes that generate pattern 
B for [A. tubingensis and A. niger). Though the pattern generated for 
A. japonicus is not identical to pattern B, it is very similar and 
cannot be distinguished clearly from it in a gel electrophoresis due 
to similar sizes and the techniques' limited resolution. A similar 



conclusion for unisereriate/biseriate species may be also deduced 
when evaluating by gel electrophoresis a profile generated after 
Rsal digestion. 

Sequence alignment with reference strain sequences retrieved 
from NCBI, together with BLAST analysis verified the RFLP 
results and showed that 3 out of the ten individual isolates were 
assigned to A. carbonarius, 3 to A. tubingensis, 2 showed high 
homology to members of the A. niger aggregate (A. niger /A. awamori) 
and 1 showed high homology to A. japonicus! A. uvarum species. 
Isolate Ai-1 was assigned to the species A. ibericus after BLAST 
analysis, showing 100% homology when aligned with A. ibericus 
AY656624.1, AY656623.1 and AY656622.1 strains. A ibericus was 
first reported by Serra et al. [61] as a new Aspergillus species 
isolated from grapes and to our knowledge this is the first report 
that this species is isolated in Greece. Furthermore, as shown from 
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Aj-1 GCCTTCGGGCGGC 
At- 5 9 GCCGCTTGTCGGC 
At- 10 GCCGCTTGTCGGC 
At- 46 GCCGCTTGTCGGC 
An- 2 GCCGCTTGTCGGC 
An- 4 GCCGCTTGTCGGC 
Ai-1 GCCGCTTGTCGGC 
Ac- 28 GCCGCTTGTCGGC 
Ac- 36 GCCGCTTGTCGGC 
Ac- 5 7 GCCGCTTGTCGGC 
*** * ♦*♦* 



CCGGGGCCTG 

C GC C G GGGGGGCGC C T 
C GC C G GGGGG GCGC C T 
C GC C G GGGGGGCGC C T 
C GC C G GGGGGG C6C C T 



C GC C G GGGGG GCGC C T 
C GC C G GGGGGGcHc T 
CGCCGGGGGGGcHcT 
CGCCGGGGGGGcHcT 
CGCCGGGGGGGcJcT 
* * * 



C 

TTGCC 
TTGCC 
TTGCC 
TGCC 
TGCC 
TGCC 
TGCC 
TGCC 
TGCC 
* 



CCCCGGGACCG-CGCC 
CCCCGGGCCCGTGCCC 
CCCCGGGCCCGTGCCC 
CCCCGGGCCCGTGCCC 
CCCCGGGCCCGTGCCC 
CCCCGGGCCCGTGCCC 
Ccic GG GC C CGCGCC C 
CclcGGGCCCGTGCCC 
CclcGGGCCCGTGCCC 
Cc|cGGGCCCGTGCCC 
+* **+* +** ** 



Hhal 



Aj-1 CGCCGGASACCCCAATGGAACACTGTCTGAAAGCGTGCAGTCTGAGTTGA 

At- 59 GCCGGAGACCCCAACACGAACACTGTCTGAAAGCGTGCAGTCTGAGTTGA 

At- 10 GCCGGAGACCCCAACACGAACACTGTCTGAAAGCGTGCAGTCTGAGTTGA 

At- 46 GCCGGAGACCCCAACACGAACACTGTCTGAAAGCGTGCAGTCTGAGTTGA 

An- 2 GCCGGAGACCCCAACACGAACACTGTCTGAAAGCGTGCAGTCTGAGTTGA 

An- 4 GCCGGAGACCCCAACACGAACACTGTCTGAAAGCGTGCAGTCTGAGTTGA 

Ai-1 GCCGGAGACACCAACACGAACACTGTCTG 

Ac- 28 GCCGGAGACACCAACACGAACACTGTCTG 

Ac- 36 GCCGGAGACACCAACACGAACACTGTCTG 

Ac- 57 GCCGGAGACACCAACACGAACACTGTCTG 

* * **^*_. *************** **** ********* ** 

Hinfl 




gtct gagtM ga 
gtctgagtIga 



GTCTGAGTlGA 

gtctgagtIga 



1^ 



Figure 5. Alignment of tlie ribosomal 5.8S-ITS amplicon sequences of representative Aspergilli isolates. A part of the amplicon 
sequences (nucleotides 101 to 200) alignment is presented. Restriction sites are presented as bold-underlined, and variable nucleotides are 
highlighted. The extra Hha\ restriction site present in the Ai-1 sequence is also shown. 
doi:1 0.1 371 /journal.pone.0093923.g005 



Table 2. 5.8S-ITS RFLP DNA fragment sizes profiles. 



Isolate 


Rest. Endonuclease 








Hinfl 


Hha\ 


Rsa\ 


Ac-28, Ac-36, Ac-57 (A. carbonarius) 


288, 191, 110, 8 


330, 178, 89 


521, 76 


At-10, At-46, At-59 {A. tubingensis) 


302, 289, 8 


207, 1 78, 1 24, 90 


599 


An-2, An-4 {A. niger aggregate) 


302, 289, 8 


207, 1 78, 1 24, 90 


523, 76 


Ai-1 {A. ibericus) 


289, 192, 110, 8 


186, 178, 145, 90 


523, 76 


Aj-1 {A. japonicus) 


293, 274, 8 


185, 178, 137, 75 


497, 78 


A. carbor)arm AF459734 


288, 191, 110, 8 


330, 1 78, 89 


521, 76 


A. ibericus AY656624 


289, 192, 110, 8 


186, 178, 145, 90 


523, 76 


A. brasiiiensis AJ280010 


290, 192, 110, 8 


207, 178, 124, 91 


524, 76 


A. niger AJ223852 


301, 290, 8 


207, 178, 123, 91 


524, 75 


A. tubingensis FJ629367 


302, 289, 8 


207, 1 78, 1 24, 90 


599 


A. awamori AM087614 


302, 289, 8 


207, 1 78, 1 24, 90 


523, 76 


A. foetidus AY373850 


302, 286, 8 


207 1 78, 1 24, 87 


520, 76 


A. aculeatus AJ279997 


273, 183, 110, 8 


185, 177, 137, 75 


496, 78 


A. japonicus FJ629335 


293, 274, 8 


185, 178, 137, 75 


497, 78 


A. uvarum AM745757 


293, 273, 8 


185, 177, 137, 75 


496, 78 



Individual isolates and reference strain isolates from grapes are clustered according to their RFLP profile deduced after sequencing of the 5.8S-ITS amplicon. Fragment 
sizes generated respect to the restriction endonuclease site position in the sequence, are given in bp. 
doi:1 0.1 371/journal.pone.0093923.t002 
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A. japonicus FJ629335 
Aj-1 

A. uvarum AM745757 
A. violaceofuscus FR733805 
A. aculeatus AJ279997 



A. fijiensis HESUOTi 
'— A. aculeatinus KFS^505^ 



A. saccharolyticus HM853552 



A. homomorphus FJ629334 



A. carbonarius AF459734 
A. lacticoffeatus FJ629336 
A. sclerotioniger FJ629353 
Ac-57 
Ac-28 
Ac-36 

A. ibericus AY656624 
J — A. sclerotiicarbonarius E\J^592^6 
Ai-1 

- A. brasiliensis AJ28001 0 
An-4 

- A. niger AJ223852 
A. tubingensis FJ629352 
A. acicfus JF450778 
A. costaricaensis FJ629326 
A. coreanus F J491 684 
At-59 

A. eucalypticola KC796400 
A. piperis FJ629352 
A. vadensis FJ629368 
At-10 

A. neoniger KC796401 
At-46 

A. foetidus AY373850 
A. awamori AN10&76M 
An-2 

A. ellipticus AY656631 
— A. heteromorphus AJISOOI 3 



Figure 6. Neighbour-Joining phylogenetic tree based on divergences of ribosomal 5.8S-ITS sequences. The alignment of ribosomal 
5.8S-ITS sequences from 10 isolates and 10 reference strains of Aspergillus spp. was performed using the CLUSTALQ program. Nucleotide divergences 
were estimated according to the Jukes-Cantor model. Node numbers represent the frequency (proportion) with which a cluster appears in 1000 
bootstrap runs. 

doi:10.1371/journal.pone.0093923.g006 



a comparLson of the A. ibericus 5.8S-ITS DNA sequences available 
at NCBI (data not shown), the fact that an extra Hhal site is 
present in all of these sequences makes it useful for characteriza- 
tion of ^. ibericus by PCR-RFLP using this restriction endonucle- 
ase. The phylogenetic analysis performed using reference species' 
sequences as outgroups supports the sequence data and is in 



accordance with data presented by other authors using the ITS 
region for phylogenetic analysis of black aspergilli [37]. 

According to various authors [20,34,37,38,58,62,63,] black 
aspergilli that have been isolated from grapes fall into the biseriate 
A. carbonarius, A. tubingensis, A. niger, A. aivamori, A. foetidus, A. 
brasiliensis and A. ibericus, and the uniseriate A. japonicus, A. uvarum 
and A. aculeatus species. Characterization of these isolates by 
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molecular means has been based on different approximations 
using RFLP [25,31,33,39,64], AFLP [30,32,36,65], or multilocus 
analysis [62,66]. RFLP and AFLP have been used successfully for 
characterization of the main ochratoxigenic aspergUli isolated 
from grapes [25,30,31,32,33,36,39,64,65]. Multilocus analysis 
based on calmodulin, beta-tubulin and ITS sequences has proved 
a useful tool for extended characterization and comparative study 
of isolates from a great variety of sources such as grapes, coffee, 
soil, surfaces, or unknown origin [37,62,66]. In this work we have 
followed a similar approximation using ITS PGR-RFLP in order 
to characterize black aspergUli isolates from grapes. We have been 
able to characterize the isolates belonging to the A. carbonarius, A. 
tubingensis and A. ibericus species, since the RFLP patterns (Table 2) 
can clearly differentiate these species from all other black aspergUli 
species isolated from grapes. 

The Aj-1 isolate could not be characterized only by molecular 
means since there is no difference between the RFLP patterns and 
ITS region DNA sequence for either species. According to the 
morphological observations the Aj-1 isolate was characterized as 
A. japonicus since its conidial size (4.3 ^m) faUs in the range of ^. 
japonicus (4—5 |im) and not in the range of^. uvarum (3-4 )Xm) [34]. 
PCR-RFLP when used in combination with other means such as 
use of morphological criteria is a reliable approximation for black 
aspergUli characterization, where, for example, differentiation of 
aspergUli as biseriate/uniseriate can initially short the strains to 
two major groups. Isolates An-1 to An-4 have been characterized 
as A. niger aggregate members due to the limitation of the ITS 
PCR-RFLP to differentiate between A. niger, A. awamori and A. 
foetidus species, thus resulting to a partial molecular characteriza- 
tion of these four isolates. 

With regard to the field study, Greece is a country with 
extended coast, although, areas of grape planting are character- 
ized by different climatic characteristics. Vineyards localisation 
can vary from very low altitude, near seaside, either arid or wet, to 
very high altitudes and wet climate regions. Referring to the 
present study, influence of the geographic localization of the 
vineyards on the incidence of black aspergUli in grapes was 
significant. The vineyards selected for this study are representative 
of these different climatic profUes existing in the country. Most 
studies that correlated the geo-climatic conditions with the 
presence and distribution of black aspergUli reported a trend of 
higher isolation incidence in regions with higher mean tempera- 
tures, but more wet areas [19,53,66,67]. These reports are partly 
in accordance with the present study, since the geographical 
localisation of sampled vineyards influenced black aspergUli 
presence, but with the higher incidences found in arid areas. In 
accordance with our results are the findings of Serra et al. [54,56] 
for Spain where higher incidences for black aspergUli were 
reported from hot and dry regions. In addition, most of the 
ochratoxigenic isolates came from regions with arid climate, Crete 
and Attica prefecture, but this could be also correlated with the 
low altitude of these vineyards in contrast with the Peloponnese 
and Macedonia where altitudes of grape cultivation are higher. 
Moreover, the incidence of black aspergiUi and the A. carbonarius 
isolation were much lower in Macedonia where temperatures are 
lower and the climate is wet compared with the other regions. In 
contrast, Chiotta et al. [65] have reported higher incidences of ^. 
carbonarius from vineyards of higher altitudes. In addition, Perrone 
et al. [20] and Visconti et al. [24] denoted the strong correlation of 
relative humidity and rainfall with A. carbonarius presence and 
OTA accumulation in grapes. FinaUy, many works denote that, 
apart from the climate of a wider region, very important role have 
also the different practices used in grape cultivation and the 
specific microclimate of the vineyard [68,69,70]. 
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Figure 7. Correlation of HPLC and ELISA values for OTA determination. Linear correlation and for the two methods assayed for OTA 

quantification (ELISA & HPLC). 
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A. niger aggregate and A. carbomrius were equally distributed in 
the present work, while only one uniseriate species was isolated 
and identified as A. japonicus. Similar distributions of A. carbonariiis 
within the A. section Migri group has been also reported by Bau et 
al. [69] for grapes from the Mediterranean Spanish coast and 
Tjamos et al. [51] for Rhodes island in Greece. These results are 
in accordance with our findings, since Crete, which has an insular 
Mediterranean climate, similar to those of Rhode and the 
Mediterranean Spanish coast, presented the higher incidence of 
A. carbonarius among black aspergiUi. In contrast. El Khoury et al. 
[71] for Lebanon and Chiotta et al. [65,66] for Argentina have 
reported relatively lower participation of ^. carbonarius within black 
aspergUli, while BattUani et al. [19] reported very high occurrence 
(20%) of uniseriate species in Italian vineyards during grape 
ripening. 

Several studies from countries all over the world during the last 
years have identified A. section Migri group as the responsible fungi 
for OTA contamination of grapes and derived products. 
Specifically, A. carbonarius is considered as the main OTA 
contamination source, and to a lesser extend species of A. niger 
aggregate, mainly A. niger and A. tubingensis. This conclusion derives 
from the higher ocliratoxigenic ability of ^. carbonarius in contrast 
to A. niger aggregate isolates that, even when they are more 
abundant within the A. section Mgri group, they are either totally 
incapable of toxin production or their toxigenic ability is scarce 
and weak [19,53,72]. In our study, all A. carbonarius isolates, with 
one exception, were capable for OTA production with most of 
them having a high ochratoxigenic potential, while from A. niger 
aggregate only one isolate identified as Aspergillus tubingensis was 
capable to produce OTA in vitro (Table 3). Similar results have 
been reported for other Mediterranean countries [19,54,56,69,73], 
although there are reports of higher ochratoxigenic incidences 
among A. niger aggregate and uniseriate spp. [39]. Ranges of 



ochratoxigenic potential for A. carbonarius isolates varied from 0 to 
34 ppm (/(g OTA g ' CYA), with most of them producing 
between 1 and 1000 ppb (Table 3). Our results are also in 
agreement with studies contacted in Greece [51,68], although 
Greek isolates with higher ochratoxigenic potential have been 
reported [31]. Particularly, Tjamos et al. [51,68] in their studies 
have also shown that most A. carbonarius isolates were strong OTA 
producers with capabilities greater than 25 ppb, while only 6% of 
the A. niger aggregate presented high OTA potential. 

Conclusions 

From our results there is a clear correlation of A. section Mgri 
group with the arid climate of Crete and Attica, whereas for colder 
and wet areas its presence is scarcer. Moreover, higher incidence 
of ochratoxigenic isolates and OTA production was recorded for 
the A. carbonarius isolates, confirming in this way that A. carbonarius 
is the main species responsible for OTA contamination in Greek 
grapes. Regarding PCR-RFLP our data indicate that A. ihericus can 
be differentiated from other aspergilli using the Hhal restriction 
endonuclease. 
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